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Abstract 

The aim of this paper is to summarize the general status of our under- 
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1 Introduction 



In this paper we present a summary of the workshop on small-x parton dynamics 
held in Lund in the beginning of March 2001. During two days we went through 
a number of theoretical and phenomenological aspects of small-x physics in short 
talks and long discussions. Here we will present the main points of these discussions 
and try to summarize the general status of the work in this field. 

For almost thirty years, QCD has been the theory of strong interactions. Although 
it has been very successful, there are still a number of problems which have not 
been solved. Most of these have to do with the transition between the perturbative 
and non-perturbative description of the theory. Although perturbative techniques 
work surprisingly well down to very small scales where the running coupling starts 
to become large, in the end what is observed are hadrons, the transition to which 
is still not on firm theoretical grounds. At very high energies another problem 
arises. Even at high scales where the running coupling is small the phase space 
for additional emissions increases rapidly and makes the perturbative expansion ill- 
behaved. The solution to this problem is to resum the leading logarithmic behavior 
of the cross section to all orders, thus rearranging the perturbative expansion into 
a more rapidly converging series. 

The DGLAP[ H, ^ ^, ^ evolution is the most familiar resummation strategy. Given 
that a cross section involving incoming hadrons is dominated by diagrams where 
successive emissions are strongly ordered in virtuality, the resulting large loga- 
rithms of ratios of subsequent virtualities can be resummed. The cross section 
can then be rewritten in terms of a pro cess- dependent hard matrix element con- 
voluted with universal parton density functions, the scaling violations of which 
are described by the DGLAP evolution. This is called coUinear factorization. Be- 
cause of the strong ordering of virtualities, the virtuality of the parton entering 
the hard scattering matrix element can be neglected (treated coUinear with the 
incoming hadron) compared to the large scale . This approach has been very 
successful in describing the bulk of experimental measurements at lepton-hadron 
and hadron-hadron colliders. 

With HERA, a new kinematic regime has opened up where the very small x parts 
of the proton parton distributions are being probed. The hard scale, Q^, is not 
very high in such events and it was expected that the DGLAP evolution should 
break down. To some surprise, the DGLAP evolution has been quite successful in 
describing the strong rise of the cross section with decreasing x. For some non- 
inclusive observables there are, however, clear discrepancies as summarized below 
in table |l]. 

At asymptotically large energies, it is believed that the theoretically correct de- 
scription is given by the BFKL [ §, §|, evolution. Here, each emitted gluon is 
assumed to take a large fraction of the energy of the propagating gluon, {\ — z) for 
2; — > 0, and large logarithms of 1/z are summed up to all orders. Although the rise 
of F2 with decreasing x as measured at HERA can be described with the DGLAP 
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Table 1: Summary of the ability of the collinear and k j^- factorization approaches to 
reproduce the current measurements of some observables: OK means a satisfactory 
description; 1/2 means a not perfect but also a not too bad description, or in part 
of the phase space an acceptable description; OK? means satisfactory description 
if a heavy quark excitation component is added in leading order; NO means that 
the description is bad; and ? means that no thorough comparison has been made. 



evolution, a strong power-like rise was predicted by BFKL. Just as for DGLAP, 
it is possible to factorize an observable into a convolution of process-dependent 
hard matrix elements with universal parton distributions. But as the virtuality 
and transverse momentum of the propagating gluon are no longer ordered, the 
matrix elements have to be taken off-shell and the convolution is also over trans- 
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verse momentum with unintegrated parton distributions. We therefore talk about 
/c_L-factorization [ |35|, ^ or the semihard approach [ p8| . 



Recently, the next-to-leading logarithmic (NLL) corrections to the BFKL equation 
were calculated and found to be huge [ |39|, This is related to the fact that 
at any finite energy, the cross section will also get contributions from emissions of 
gluons which take only a small fraction of the energy of the propagating gluon. 



The CCFM [ ^ ^ evolution equation resums also large logarithms of 

1/(1 — z) in addition to the 1/z ones. Furthermore it introduces angular ordering 
of emissions to correctly treat gluon coherence effects. In the limit of asymptotic 



energies, it is almost equivalent to BFKL [ J5, 56, J7|, but also similar to the 



DGLAP evolution for large x and high Q^. The cross section is still /c^-factorized 
into an off-shell matrix element convoluted with an unintegrated parton density, 
which now also contains a dependence on the maximum angle allowed in emissions. 

An advantage of the CCFM evolution, compared to the BFKL evolution, is that 
it is fairly well suited for implementation into an event generator program, which 
makes quantitative comparison with data feasible also for non-inclusive observables. 
There exist today three such generators [ |g, |0|, |1], |2|, |3|, |4|, m and 



they are all being maintained or/and developed by people from the departments 
of physics and of theoretical physics at Lund University. 

Since 1998 there has been an on-going project in Lund, supported by the Swedish 
Royal Academy of Science, trying to get a better understanding of the differences 
between the different generators and to compare them to measured data. This 
project is what led up to the meeting in Lund in early March 2001, where a number 
of experts in the field were invited to give short presentations and to discuss the 
current issues in small-x physics in general and /c^-factorization in particular. 

In this article we will try to summarize these discussions and give a general status 
report of this field of research. We also suggest the formation of an informal 
collaboration of researchers in the field, to facilitate a coherent effort to solve some 
of the current problems in small-x parton dynamics. 

The outline of this article is as follows. First we give a general introduction to 
fc^-factorization in section |^. Then, in section ^ we discuss the off-shell matrix 
elements, both at leading order and the prospects of going to next-to-leading order. 
In section ^, we discuss the unintegrated parton distributions and how they are 
evolved. Here we also try to quantify the transition between DGLAP and BFKL. 
We present a number of parameterizations of unintegrated parton distributions 
and make a few comparisons. Then we describe the next-to-leading logarithmic 
correction to the evolution. In section |^ we describe the available event generators 
for small-x evolution. Finally in section ^ we present our conclusions and discuss 
the forming of an informal collaboration to better organize the future investigations 
of small-x phenomenology. 
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Q = {-xP+,Q-,6) 




P=(P+, 0,0) 

Figure 1: Schematic picture of a typical unitarity diagram for deep inelastic scatter- 
ing. An incoming proton with a large positive light-cone momentum, , is being 
probed by a photon with a large virtuality and a large negative light-cone momen- 
tum. The photon scatters on a parton from the proton with space-like momentum 
k. 



2 The fc^-fcictorization approach 

The calculation of inclusive quantities, like the structure function F2{x,Q^) at 
HERA, performed in NLO QCD is in perfect agreement with the measurements. 
The NLO approach, although phenomenologically successful for F2{x,Q'^), is not 
fully satisfactory from a theoretical viewpoint because, in the words of Catani, "the 
truncation of the splitting functions at a fixed perturbative order is equivalent to 
assuming that the dominant dynamical mechanism leading to scaling violations is 
the evolution of parton cascades with strongly ordered transverse momenta" [ 

As soon as exclusive quantities like jet or heavy quark production are investigated, 
the agreement between NLO coefficient functions convoluted with NLO DGLAP 
parton distributions and the data is not at all satisfactory. Large so-called K- 
factors (normalization factors, for example K = ) are needed to bring the 
NLO calculations close to the data [ 0, |5^, ^ {k ~ 50 for J/ip production 
and K 2 — 4 for bottom production at the TEVATRON), indicating that a 
significant part of the cross section is still missing in the calculations. 

At small X the structure function F2{x,Q'^) is proportional to the sea quark den- 
sity, which is driven by the gluon density. The standard QCD fits determine the 
parameters of the initial parton distributions at a starting scale Qo- With the help 
of the DGLAP evolution equations these parton distributions are then evolved to 
any other scale Q^, with the splitting functions still truncated at fixed 0{as) (LO) 
or (NLO). Any physics process in the fixed order scheme is then calcu- 
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Figure 2: Kinematic variables for multi-gluon emission. The t-channel gluon mo- 
menta are given by ki and the gluons emitted in the initial state cascade have 
momenta pi. The upper angle for any emission is obtained from the quark box, as 
indicated with E. We define z±i = k±i/k±(^izpi) and qi = p±i/{l — 



lated via collinear factorization into the coefficient functions C"(|) and coUinear 
(independent of k±) parton density functions: fa{z,Q'^): 



a 



ao [-C^i-)faiz,Q' 
J z z 



At large energies (small x) the evolution of parton distributions proceeds over a 
large region in rapidity Ay ~ log(l/x) and effects of finite transverse momenta of 
the partons may become increasingly important. Cross sections can then be k± - 
factorized [ ^ |3^, ^ into an off-shell {k± dependent) partonic cross section 



a{-, kj_) and a. k± - unintegrated parton density function J-'{z, kj 



a= [ —d'^k±a( — ,k\)J^{z,k'] 
J z z 



(2) 



The unintegrated gluon density J^{z, /c^) is described by the BFKL [ ||, |], 0] 
evolution equation in the region of asymptotically large energies (small x). An 
appropriate description valid for both small and large x is given by the CCFM 
evolution equation [ ^ resulting in an unintegrated gluon density 



A{x,kj_,q'^), which is a function also of the additional scale q described below. 
Here and in the following we use the following classification scheme: xQ{x, kj_) de- 
scribes DGLAP type unintegrated gluon distributions, xJ-'{x, kj_) is used for pure 
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BFKL and xA{x,kl,f) stands for a CCFM type or any other type having two 
scales involved. 

By explicitly carrying out the k± integration in eq.(^ one can obtain a form fully 
consistent with coUinear factorization [ |5^, ^ : the coefficient functions and also 
the DGLAP splitting functions leading to fa{z, Q^) are no longer evaluated in fixed 
order perturbation theory but supplemented with the all-loop resummation of the 
Oglogl/x contribution at small x. This all-loop resummation shows up in the 
Regge form factor Ajiegge for BFKL or in the non-Sudakov form factor A^^ for 
CCFM, which will be discussed in more detail in chapter ^ 

3 OfF-shell matrix elements 



LO NLO resolved photon 




fcx = fcx = 




k± — factorisation 



Figure 3: Diagrammatic representation of LO, NLO and resolved photon processes 
in the collinear approach (top row) and compared to the k±- factorization approach. 

It is interesting to compare the basic features of the fc^-factorization approach to 
the conventional collinear approach. In /c^-factorization the partons entering the 
hard scattering matrix element are free to be off-mass shell, in contrast to the 
collinear approach which treats all incoming partons massless. The full advantage 
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of the fcx-factorization approach becomes visible, when additional hard gluon radi- 
ation to a 2 ^ 2 process like 7(7 QQ is considered. If the transverse momentum 
P±g of the additional gluon is of the order of that of the quarks, then in a con- 
ventional coUinear approach the full O(a^) matrix element for 2 ^ 3 has to be 
calculated. In fc^-factorization such processes are naturally included to leading log- 
arithmic accuracy, even if only the LO off-shell matrix element is used, since the 
k± of the incoming gluon is only restricted by kinematics, and therefore can aquire 
a virtuality similar to the ones in a complete fixed order calculation. In Fig. |^ we 
show schematically the basic ideas comparing the diagrammatic structure of the 
different factorization approaches. Not only does /c^-factorization include (at least 
some of the) NLO diagrams [ |60| it also includes diagrams of the resolved photon 
type, with the natural transition from real to virtual photons. 

However, it has to be carefully investigated, which parts of a full fixed NLO cal- 
culation are already included in the /c^-factorization approach for a given off-shell 
matrix element and which are still missing or only approximatively included. It 
should be clear from the above, that the 0{as) matrix element in fc^-factorization 
includes fully the order 0{as) matrix element of the collinear factorization ap- 
proach, but includes also higher order contributions. In addition, due to the unin- 
tegrated gluon density, also parts of the virtual corrections are properly resummed 
(Fig.|b,c). 



3.1 Order off-shell matrix elements 

Several calculations exist for the process •yg — > QQ and gg — > QQ, where the gluon 
and the photon are both allowed to be off-shell [ |3^ and Q{Q) can be a heavy 
or light quark (anti quark). 




(a) (6) (c) 



Figure 4: Schematic diagrams for k±-factorization: (a) shows the general case for 
hadroproduction of (heavy) quarks, (b) shows the one-loop correction to the Born 
diagram for photoproduction (c) shows the all-loop improved correction with the 
factorized structure function J-'{x, kjj or A{x, k'j_, (f). 
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The four- vectors of the exchanged partons ki, k2 (in Sudakov representation) are 
(Fig|a): 

fcf = ZiPf + z,P,^ + k'i, (3) 
k^ = z^Pr + Z2P^ + A;l2 (4) 

with Zj, Zi being the two components (+, — ) of the hght cone energy fraction. The 
(heavy) quark momenta are denoted by ^3, and the incoming particles (partons) 
by Pi, P2 with 

^^1,2 = ^7^(0, ±1,1), 2P,P2 = S (5) 

where indicates the vanishing two-dimensional transverse momentum vector. In 
the case of photopro duct ion or leptopro duct ion this reduces to: 

ki = Pi, = photo-production (6) 
ki = = yPi + yP2 + q±2^ = leptoproduction (7) 

In all cases the off-shell matrix elements are calculated in the high energy approx- 
imation, with Zi = Z2 = 

fcf = ziP!^ + k^^i (8) 
ki^ = Z2P^ + A;^2, (9) 

which ensures that the virtualities are given by the transverse momenta, k\ = —kj_i 
and = ~kj_2- The off-shell matrix elements involve 4-vector products not only 
with ki,k2 and the outgoing (heavy) quark momenta P3,Pa, but also with the 
momenta of the incoming particles (partons) Pi, ^2- This is a result of defining the 
(off-shell) gluon polarization tensors in terms of the gluon (ki) and the incoming 
particle vectors (Pi, 2); which is necessary due to the off-shellness. In the coUinear 



limit, this reduces to the standard polarization tensors. In [ Q it has been 
shown analytically, that the off-shell matrix elements reduce to the standard ones 
in case of vanishing transverse momenta of the incoming (exchanged) partons ki, 

k2. 

Due to the complicated structure of the off-shell matrix elements, it is also necessary 
to check the positivity of the squared matrix elements in case of incoming partons 
which are highly off shell [kl^k^ <^ 0). It has been proven analytically in [ |61 



for the case of heavy quarks with both incoming partons being off mass shell. 
We have also checked numerically, that the squared matrix elements are positive 
everywhere in the phase space, if the incoming particles (electron or proton) are 
exactly massless (Pi^ = P2 =0). As soon as finite masses (of the electron or proton) 
are included, the exact cancellation of different terms in the matrix elements is 
destroyed, and unphysical (negative) results could appear. 

From the off-shell matrix elements for a 2 ^ 2 process, it is easy to obtain the high 
energy limit of any on-shell 2(massless partons) 2+n(massless partons) process, 
with n = 1,2 (see Fig. ^jb for the case of photoproduction), by applying a correction 



8 



as given in [R^, page 180]. By doing so, the corrections coming from additional real 
gluon emission are estimated (i.e. 7(7 QQg corrections to ■yg QQ) with the 
incoming gluons now treated on-shell. In the case of hadropro duct ion one obtains 
the NNLO corrections gg QQgg- In this sense /c_|_-factorization provides an easy 
tool for estimating i^-factors in the high energy limit: 



K 



(7 tot 



a{LO or NLO) 



a(2 ^ 2) +a(2 ^ 2 + 
cr(2 2) 



(10) 



In Tab. ^ we give an overview over the different off-shell matrix elements avail- 
able. It has been checked, that the different calculations [ |3^, |62| of the pro- 
QQ give numerically the same result. Also the matrix elements of [ 



cess 'y^*^g 



35|, p6|, Bl] are identical. 



process 


comments 


reference 


ig* ^ QQ 
i*g* ^ QQ 
i*g* ^ QQ 
g*g* ^ QQ 


]M]2 
]M]2 

partially integrated 
]M]2 


[ 35, 62, 19] 
[ 35, 62] 
64] 


ig* J/i^g 

7^* ^ J/^g 

Yg* J/i^g 

g*g* J/ipg 


]M]2 

helicity amplitude 
helicity amplitude 
helicity amplitude 


[0, 3C] 

m 
m 
m 



Table 2: Table of calculations ofO{as) off-shell matrix elements for different types 
of processes. 



3.2 Next-to-Leading corrections to OLg ofF-shell matrix ele- 
ments 

So far we have been dealing with matrix elements only in 0{as), but even here 
/c_L-factorization has proven to be a powerful tool in the sense that it includes al- 



ready in lowest order a large part of the NLO corrections [ 6C] of the coUinear 
approach. However, LO calculations are still subject to uncertainties in the scale 
of the coupling and in the normalization. As the next-to-leading corrections (in- 
cluding energy-momentum conservation) to the BFKL equation are now available 
also the matrix elements need to be calculated in the next-to-leading order. 

NLO corrections to the process y*g* — ^ X contain the virtual corrections to 7* (7* 
qq (Fig. ^a) and the leading order off-shell matrix element for the process 7*5'* 
qqg (Fig. ^). The calculation of the former ones has been completed, and the 
results for the matrix elements are published in [ The real corrections (7*5'* 



qqg) have been obtained for helicity-summed squares of the matrix elements. For 
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+ 




Figure 5: Schematic diagrams for next-to-leading contributions: (a) virtual cor- 
rections, (b) real corrections, (c) the process 7* + g — > (qq) + q, (d) the process 
'y* + q ^ (qqg) + q, (e) diagrams showing the reggeization of the gluon 



the longitudinally polarized photon they are published in [ and results for the 
transversely polarized photon will be made available soon [ 170 . 



The starting point of these calculations is a study of the Regge limit of the processes 
Y + Q ~^ (qq) + Q (Fig- ib) 7* + ^ {<1<19) + Q (Fig. i.e. the contributing 
QCD diagrams have been calculated in the region s = {q + pY ^ Q^, M^, t and 
^ Aqcd, with q{p) being the four- vector of the photon (quark). In this limit, 
the exchanged gluon is not an elementary but a reggeized gluon: the two gluon 
exchange diagrams (Fig. |^) contain a term proportional to uj{t) log(s), where u;(t) 
is the gluon Regge trajectory function. This term is not part of the subprocess 
'-f*g* qq but rather belongs to the exchanged gluon. In order to find the con- 
tribution of Fig. 1^ to this scattering subprocess, one first has to subtract the 
reggeization piece. This means that, for the subprocess 'y*g* — > qq, the notation 
in NLO not only stands for 'off-shellness' but also for 'reggeized'. It also has 
important consequences for the factorization of the NLO off-shell subprocess inside 
a larger QCD diagram: the t-channel gluons connecting the different subprocesses 
(see, for example, g{ki) in Fig. ^) are not elementary but reggeized, and the QCD 
diagrams include subsets of graphs with more than two gluons in the t-channel (see 
Fig. |e) 

The results for the virtual corrections contain infrared singularities, both soft and 
virtual ones. As usual, when integrating the real corrections over the final state 
gluon, one finds the same infrared singularities but with opposite signs. So in 
the sum of virtual and real corrections one obtains an infrared finite answer. A 
peculiarity of this NLO calculation is the interplay with the LO process 7* + 
q (qq) + g + q. In the latter, the process is calculated in the leading log(s) 
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approximation (LO BFKL approximation), where the gluon has a large rapidity 
separation w.r.t. the gg-pair (i.e. it is emitted in the central region between the 
gg-pair and the lower quark line in Fig. |^). As a result of this high energy (small 
X- or Regge) factorization, expressions for the vertex g* + g* ^ g and for the 
subprocess 7*5'* — qq are obtained. When turning to the NLO corrections of 
the process 7* + g — > (qqg) + g, the calculation extends to next-to-leading log(s) 
accuracy, but it contains, as a special case, also the LO BFKL process. In order to 
avoid double counting, one has to subtract the central region contribution. Only 
after this subtraction we have a clean separation: qqg final states without or with 
a large rapidity separation between the gluon and the quark-antiquark pair. The 
former one belongs to the NLO approximation, whereas the latter one is counted 
in the leading log(s) approximation. 



The results for real corrections in [ ^ are very interesting also in the context of 



the QCD color dipole picture. It is well-known that the total LO 7*g cross section 



at high energies can be written in the form [ H, 72 



= dz dp {Q, z, p)) agg{x, z, p)^p^*g (Q, z, p) (11) 



where z and (1 — z) denote the momentum fractions (w.r.t. the photon momentum 
g) of the quark and the antiquark, respectively, p is the transverse extension of the 
quark-antiquark pair, ipqq stands for the gg Fock component of the photon wave 
function, and aqq is the color dipole cross section. This form of the scattering cross 
section is in agreement with the space time picture in the target quark rest frame: 
a long time before the photon reaches the target quark at rest, it splits into the 
quark-antiquark pair which then interacts with the target. During the interaction 
the transverse extension of the quark-antiquark pair remains frozen, i.e the initial 
quark-antiquark pair has the same p-value as the final one. Via the optical theorem 
the 7*g total cross section is related to the square of the scattering matrix element 
of the process 7* + g {qq)+q', the dipole cross section form eq.(|n|) must then be a 
consequence of the special form of the LO amplitude for the subprocess Yo* ~^ QQ- 
Since this form of the total 7*g cross section (which is preserved when the target 
quark is replaced by the target proton) looks so appealing (and also has turned 
out to be extremely useful in phenomenological applications), that it is desirable 
to investigate its validity also beyond leading order. 

When trying to generalize eq.(pJ^ to NLO, one is led to study the square of the 
real corrections 7* + g — ^ {QQg) + Q- If the color dipole picture remains correct 
(factorization in wave function and dipole cross section), this contribution naturally 
should lead to a new qqg Fock component of the photon wave function, and to a 
new interaction cross section (Tqqg, which describes the interaction of the quark- 



ant iquark-gluon system with the target quark. In [ |6^, |70[ it is shown that this 
is indeed the case. The form of the new photon wave function is rather lengthy 
(in particular for the transverse photon), and a more detailed investigation is still 
needed. Nevertheless, one feels tempted to conclude that the color dipole form 
eq.(pJ]) is the beginning of a general color multipole expansion, where the different 
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Fock components of the photon wave function describe the spatial distribution 
of color charge inside the photon. However, before this conclusion can really be 
drawn, it remains to be shown that also the virtual corrections to the scattering 
amplitude of 7* + g — ^ (qq) + g fit into the form eq.(|Tl|). These corrections should 
lead to NLO corrections of the photon wave function or the dipole cross section. 
They also may slightly 'melt' the transverse extension of the quark- ant iquark pair 
during the interaction with the target. 

In summary, the results in [ ^ |70| provide the starting point for discussing 



exclusive final states in the /c^-factorization scheme in NLO accuracy. However, be- 
fore these formulae can be used in a numerical analysis, virtual and real corrections 
have to be put together, and IR finite cross sections have to be formulated. The 
NLO corrections to the off-shell matrix elements described in this subsection also 
represent the main ingredients to the NLO photon impact factor. Its importance 



will be discussed in section 4.6.1 



3.3 Gauge-invariance 

The off-shell matrix elements and the cross section taken in lowest order as are 
gauge invariant, as argued in [ However, when extended to next order in 
perturbation theory, as depicted in Fig. ^ problems will occur and the off-shell 
matrix elements and unintegrated parton distributions are no longer necessarily 
gauge invariant. The following critique indicates that their definition needs further 
specification and that further work is needed to properly justify the formalism [ 

m- 

Basically, parton distributions are defined as a hadron expectation value of a quark 
and antiquark field (or two gluon fields): 

{p\il;{y)r4j{0)\p), (12) 

with an appropriate Fourier transformation on the space-time argument y^, some 
appropriate Dirac matrix F and some appropriate normalization. This definition 
is not gauge invariant, so one must either specify the gauge or one modifies the 
definition to make it gauge invariant: 

{p\ij{y)TPe-'3^o'^y'"^"^y'^'"iP{0)\p). (13) 

Here, ta are generating matrices for the SU (3) color group of QCD, and the symbol 
P denotes that the gluon fields are laid out in their order on the path in the integral 
from to y. 

The big question is which path is to be used. For the integrated parton distribu- 
tions, 1/ is at a light-like separation from 0: normally y~ 0, y~^ = y± = 0, and 
then taking the path along the light-like straight line joining and y = (0, y~ , 0±) 
is natural and correct. But for unintegrated distributions, this is not so simple; one 
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has a choice of paths. The choice is not arbitrary but is determined by the deriva- 
tion: i.e., the definition is whatever is appropriate to make a correct factorization 
theorem. 

This can be seen from the analysis of gluon emission that we have aheady de- 
scribed. This analysis is only applicable if ladder diagrams, as in Fig. |^, actually 
dominate. In fact, in a general gauge, non-ladder diagrams are as important as lad- 
der diagrams. This is the case, for example, in the Feynman gauge. The standard 
leading-logarithm analysis suggests, as is commonly asserted, that the appropriate 
gauge is a light-cone gauge n-A = 0, where n is a light-like vector in a suitable direc- 
tion, for then non-ladder contributions are suppressed in the LLA analysis. From 
this one would conclude that the off-shell matrix elements and the unintegrated 
parton distributions are defined as the appropriate field theory Green functions in 
light-cone gauge. 

Closer examination shows that there must be problems since the unintegrated par- 
ton distributions are divergent. This was shown quite generally by Collins and 
Soper [ ^ |75[ as part of their derivation of factorization for transverse-momentum 



distributions. They found they needed to derive factorization and to define un- 
integrated parton distributions in a non-light-like axial gauge, i.e., with 7^ 0.[| 
They derived an evolution equation for the dependence on the gauge-fixing vector; 
the evolution is quite important and contains doubly-logarithmic Sudakov effects. 
In the limit — > 0, the parton distributions become so singular that they are 
not defined. The divergences are associated with the emission of gluons of arbi- 
trarily negative rapidity with respect to the parent hadron; commonly these are 
called soft gluons, and a non-light-like gauge-fixing vector provides a cutoff on the 
divergences. It is normally said that soft gluons cancel in QCD cross sections. 
However, the standard cancellation requires an integral over all parton transverse 
momentum, which is of course not applicable in an unintegrated parton density. 

The divergences can readily be seen in one-loop graphs, and are a generalization of 
divergences known to occur in light-cone gauge. Modification of the ie prescription 
of the 1/k-n singularity of the gluon propagator, as is commonly advocated, is not 
sufficient, since this does not remove the divergence in the emission of rea/ gluons, 
for which the singularity is an endpoint singularity in k ■ n. 

Clearly, in all the derivations of the BFKL and related equations, there must be 
an appropriate cutoff. Unfortunately, this is not normally made very explicit, even 
though an explicit specification of the cutoff is vitally important if the formalism 
is to make sense beyond LO. One can expect, that such a cutoff is a cutoff in 
angle or rapidity, and it should be related to angular ordering, supporting the 
intuitive approximate picture. However, the presence of the cutoff implies that 
there is an extra parameter in the parton distributions, whose variation should be 



understood. Balitsky [ ^ has worked in this direction, although his formulation 



does not appear to have developed far enough to address the issues in the present 



^ An equivalent and probably better definition can be made with suitable path-ordered expo- 
nential in the operator. 
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document. 

A proper derivation will also result in an explicit definition of a reggeized gluon. 
Such a definition is not readily extracted from the original BFKL publications. At 
most these publications provide a definition as a property of a solution of their 
equation. No explicit definition is given whereby quantities involving reggeized 
gluons are matrix elements of some operator or other. 

4 Unintegrated parton distributions 
4.1 Introduction 

The conventional parton distributions describe the density of partons carrying a 
certain longitudinal momentum fraction inside the proton. These distributions are 
integrated over the transverse momenta of the partons up to the factorization scale 
/X. However, in order to describe some exclusive processes it becomes necessary to 
consider the transverse momenta of the partons and thus use so called unintegrated 
gluon distributions A{x, kj_, fi^). Unintegrated parton distributions account for the 
resummation of a variety of logarithmically enhanced terms, such as 

(a, log(/iVA2))" , («. log(/iVA2) log(l/x))" , (a, log(l/x))" and (a, \og\fxykl)) 

(14) 

The unintegrated parton distributions describe the probability of finding a parton 
carrying a longitudinal momentum fraction x and a transverse momentum k± at 
the factorization scale /i. The unintegrated gluon density xA{x, kj_, fj,"^) can be 
related to the integrated one xg{x, /i^) by: 

X5f(x, /i^)^ / dk'{xA{x,k']_, ^'^) (15) 
Jo 

The ~ sign in the above equation indicates, that there is no strict equality between 
unintegrated and integrated parton distributions, as neither are observables. 

In unintegrated parton distributions, the contribution from the large logarithms 
in X and /i^, the terms specified above in expression (p!4D, can be correctly dis- 
entangled. However, the unintegrated parton distributions are defined only as a 
functions of three variables x, k'j_, /i^. The 4- vector k of the propagator gluon in 
initial state cascade is given by: 

k = {k+,k-,k±) 
with k~^ = x~^P~^, k~ = x~ P~ The virtuality k"^ is: 

1.2 _ L+L- _ l2 _ + - p- p+ _ 1.2 

rV th th J_ Uj Jy 1 1 J_ * 

In the region of strongly ordered k~^, the typical values of k~ are small enough 
that it can be neglected in the hard scattering factor. For example, the virtuality 
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is dominated by the transverse momentum only: 



1^' ^ -A; 



However, the k~ integral is still present, and in fact, part of the definition of 
an unintegrated parton density is, that the relevant matrix element of a gluonic 
operator is integrated over all k~ . This is the reason, why only x = is kept in 
the unintegrated gluon density. 

The all-loop resummation of the log 1/x contribution at small x leads to Reggeiza- 
tion of the gluon vertex, giving rise to a significant form factor. The Regge form 
factor Apj^egge, often also called non-Sudakov form factor Ans, regularizes the 1/z 
divergence in the splitting function. 



{Z+} OC 



(16) 



The ARegge of BFKL is given by [ 0: 



with as 



^ Regge I 



kl) = exp (-as f'^ f ^Q{kl - q'Deiq'l - 
\ J z J qj_ 



(17) 



3a s 



and go being a lower cutoff. This form factor can be ex- 



panded by a power series and then symbolically represented as: 





as(^±)^ 



tts log (z) log 



^as log (z) log 



where the small x resummation of the virtual corrections becomes obvious. Such a 
diagrammatic representation is of course gauge dependent. It should be noted that 
the Regge form factor in this way looks like the Sudakov form factor [ |7^, |7S[ used 
to regularize the 1/(1 — z) pole in standard DGLAP evolution, but the resummed 
diagrams involve small z rather than large ones. This similarity is used in the 
derivation of the Linked Dipole Chain model below. 

In interpreting the results quoted in this section, the reader should bear in mind 
the caveats explained in Sec. 3.3, that as yet, no fully explicit definition of the 
parton distributions has been given. 



4.2 The CCFM evolution 



According to the CCFM [ ^ ^| evolution equation the emission of gluons 
during the initial cascade is only allowed in an angular-ordered region of phase 
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space. The maximum allowed angle H is defined by the hard scattering quark box, 
producing the (heavy) quark pair. In terms of Sudakov variables the quark pair 
momentum is written as: 

p,+P5 = T(Pi + SP2) + g± (18) 

where T (TS) is the positive (negative) light-cone momentum fraction of the quark 
pair, Pi 2 are the four-vectors of incoming beam particles, respectively and Qt is 
the vectorial sum of the transverse momenta of the quark pair in the laboratory 
frame. Similarly, the momenta pi of the gluons emitted during the initial state 
cascade are given by (here treated massless): 

2 

Pi = V,{Pl + i,P2)+P^^, ^i = ^, (19) 

svf 

with Vi = {1 — Zi)xi-i, Xi = ZiXi-i and s = (Pi -|- P2)^ being the squared center of 
mass energy. The variable C,i is connected to the angle of the emitted gluon with 
respect to the incoming proton and Xi and Vi are the momentum fractions of the 
exchanged and emitted gluons, while Zi is the ratio of the energy fractions in the 
branching {i — 1) ^ i and p±i is the transverse momentum of the emitted gluon i. 

The angular-ordered region is then specified by (Fig. pj): 

^0 < ei < ■ ■ ■ < a < S (20) 

which becomes: 

Zi-iqi-i < qi (21) 
where the rescaled transverse momenta qi of the emitted gluons is defined by: 

qi = Xi-i^/s^i = (22) 

1 Zi 

It is interesting to note, that the angular ordering constraint, as given by eq. (pTD , 
reduces to ordering in transverse momenta px_ for large whereas for z — > 0, the 
transverse momenta are free to perform a so-called random walk. 

The CCFM evolution equation with respect to the scale can be written in a 
differential form 



_2 0? xA{x,k\,q^) r d(pP{z,{q/zf,k\) , 



Q 



where Ai^x.k^.cf) is the unintegrated gluon density, depending on x, and the 
evolution variable = (f'. The splitting variable is z = x/x' and kj_ = (1 — 
z)/zq + where the vector q is at an azimuthal angle (p. The Sudakov form 
factor As is given by: 

9' dq^ /-i-Qo/g as (q^(l 



A.(f,Ql)^^,[^-j-^l -^ -7:/" I (24) 
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with as = = For inclusive quantities at leading- logarithmic order the Su- 
dakov form factor cancels against the 1/(1 — z) coUinear singularity of the splitting 
function. 

The splitting function P for branching i is given by: 

P.U. ll kid = ""'fi'^-'"' + ^A„.(.., 4.) (25) 
where the non-Sudakov form factor A^^ is defined as: 

\ogKs{z^,qlkl,) = -as f ^ f %Q{k^^ - q)Q{q - z' q,) (26) 

The upper limit of the z' integral is constrained by the functions in eg . (1261) by: 
Zi < z' < min(l, fc^j/gj), which results in the following form of the non-Sudakov 
form factor [ |8^ : 

log A„, = -as(fcL) log (-) log (27) 

\Zi/ yzoZiQi ) 

where 

r 1 if k^ijqi > 1 

Zq = 1 k±i/p±i a Zi < k±i/qi < 1 
[ Zi if kx.il qi < Zi 

The unintegrated gluon density A{x, k'j_,q'^) is a function also of the angular vari- 
able q, ultimately limited by an angle, q^ = x^.^Ss, defined by the hard interaction, 
and the two scales k'j_, (f in A{x^ k\, (f) should not come as a surprise. As we are 
concentrating on an evolution, which is not ordered in transverse momenta, it is 
natural to have more than one scale in the problem. The typical example of such 
a two scale evolution process is 7*7* — > hadrons scattering where the photons are 
highly virtual, but the virtualities still being much smaller than the total energy, 
s ^ Aqcd- Another example is DIS with a forward high-p^ jet, where 

and p\ provide the scales. In the DGLAP approximation the evolution is per- 
formed between a small and a large scale - a large scale probe resolves a target at 
a smaller scale. It is obvious, that this evolution is not appropriate for the case of 
two similar scales and instead an evolution in rapidity, or angle, is needed. 

In CCFM the scale q (coming from the maximum angle) can be related to the 
evolution scale in the coUinear parton distributions. This becomes obvious since 

q^ ~ x^Hs = yxgs = s + Q\ (28) 

The last expression is derived by using + ~ XgP2 + yPi + Qx_, 'B.'r^yjxg and 
s = yXgS — Q\. This can be compared to a possible choice of the renormalization 
and factorization scale //^ in the coUinear approach with /i^ = Q\ + 4 ■ and the 
similarity between yU and q becomes obvious. 
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4.3 LDC and the transition between BFKL and DGLAP 



The Linked Dipole Chain (LDC) model [ 50, ^ is a reformulation of CCFM 
which makes the evolution explicitly left-right symmetric, which will be discussed 



in section in more detail. LDC relies on the observation that the non-Sudakov 
form factor in eg. (p6D , despite its name, can be interpreted as a kind of Sudakov 
form factor giving the no-emission probability in the region of integration. To 
do this an additional constraint on the initial-state radiation is added requiring 
the transverse momentum of the emitted gluon to be above the smaller of the 
transverse momenta of the connecting propagating gluons: 

pL>min(fcL_i,fcL). (29) 

Emissions failing this cut will instead be treated as final-state emissions and need to 
be resummed in order not to change the total cross section. In the limit of small z 



and imposing the kinematic constraint of eq. ( |64|) (see section 4^ for further details) 
gives a factor which if it is multiplied with each emission, completely cancels the 
non-Sudakov and thus, in a sense de-reggeizes the gluon. 

One may ask when it is appropriate to use collinear factorization (DGLAP) and 
when we must account for effects of BFKL and /c^-factorization. Clearly, when k± 
is large and 1/x is limited we are in the DGLAP regime, and when 1/x is large 
and k_i_ is limited we are in the BFKL regime. An essential question is then: What 
is large? Where is the boundary between the regimes, and what is the behavior 
in the transition region? Due to the relative simplicity of the LDC model, which 
interpolates smoothly between the regimes of large and small k±, it is possible 
to give an intuitive picture of the transition. In the following x is always kept 
small, and leading terms in log 1/x are studied. Thus the limit for large k± does 
not really correspond to the DGLAP regime, but more correctly to the double-log 
approximation. 

First we will discuss the case of a fixed coupling. For large values of k± (in the 
"DGLAP region") the unintegrated structure function !F{x,k^) is dominated by 
contributions from chains of gluon propagators which are ordered in k± . The result 
is a product of factors «s ■ ^ ■ ^ [ 0, |, i 

N ^ 

T{x, k\) ~ XI / II'^s dli6{li+i - U) dKiQ{Ki+i - Ki) 

N 

where as = — / = log(l/x) and k = log(A;^/A^) (30) 

TT 

Integration over Ki with the 6-functions gives the phase space for ordered values 
Ki. The result is k^ /N\. The integrations over li give a similar result, and thus we 
obtain the well-known double-log result 

Hx. fcl) ~ E ■ ^ ■ Tn = /o(2v^) (31) 
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In the BFKL region also non-ordered chains contribute, and the result is a power- 
like increase ~ for small x- values [ ^, H, . 

The CCFM [ |lH2], ||, H and the LDC model [ |OHT], |T1 interpolate between 



the DGLAP and BFKL regimes. In the LDC model the possibility to "go down" 
in /c_|_, from k' to a smaller value k, is suppressed by a factor exp(K — k'). The 
effective allowed distance for downward steps is therefore limited to about one unit 
in K. If this quantity is called 5, the result is, that the phase space factor ^ is 

replaced by ^^^^P ■ 
Thus we obtain: 

Hx, - L (32) 

When hi is very large, this expression approaches eq.(^). When n is small we 
find, using Sterling's formula, that /c^.) ~ Hn{.^s15^)^ /N\ = exp(A/) = 
with A = as^e. For 6 = 1 this gives A = eas = 2.72as, which is very close to the 
leading order result for the BFKL equation, A = 4 log 2 as = 2.77as- Thus eq. ( P^ ) 
interpolates smoothly between the DGLAP and BFKL regimes. It is also possible 
to see that the transition occurs for a fixed ratio between k, and /, n/l ^ eas- 

For a running coupling a factor as c^^^ in eq.(^) is replaced by aodn/hz = a^du, 
where ao = as/^, u = log^K/no) and kq is a cutoff. In the large k± region (the 
"DGLAP region") the result is therefore similar to eq. (PTI), but with as^ replaced 



by aou. For small x-values we note that the allowed effective size of downward 
steps, which is still one unit in k, is a sizeable interval in u for small n, but a very 
small interval in u for larger k. 

This is the reason for an earlier experience [0, which showed that a typical chain 
contains two parts. In the first part the A;_|_-values are relatively small, and it is 
therefore easy to go up and down in k±, and non-ordered fc^-values are important. 
The second part is an ordered (DGLAP-type) chain, where k± increases towards 
its final value. 

Let us study a chain with links, out of which the first N — k links correspond to 
the first part with small non-ordered k± values, and the remaining k links belong 
to the second part with increasing k±. Assume that the effective phase space for 
each Ui in the soft part is given by a quantity A. Then the total weight for this 
part becomes A^~^. For the k links in the second, ordered, part the phase space 
becomes as in the DGLAP case /k\. Thus the total result is {G = kJF) 

G'-l.^^l.^A = L ^, L — • (33) 

This simple model also interpolates smoothly between the DGLAP and BFKL 
regions. For large w-values the last term in the sum over k dominates, which gives 
the "DGLAP" result 

G ~ ^(ao/M)^/(A^!)^ = Io{2^/aJu) ^ (levr^ao/w)"^/^ ■ exp(2y^). (34) 



19 



^ 4 
3 
2 
1 




20 



Full model 
BFKL limit 
DGLAP limit 

600 



7*10^ 







0.5 



1 



1.5 



(GeV^) 



Figure 6: Results for a running coupling for XI = 3, which for A = 0.3 means 
X = 5 ■ 10~^. (Note that this x-value corresponds only to the chain between the 
"observed" low energy gluon and the parent soft gluon.) The solid line corresponds 
to the model in eq. ^3^), the dashed line to the BFKL limit, eq. 1^^), and the dotted 
line to the DGLAP limit in eq. ^3^). The abscissa is the variable w = aou/X^l 
defined so that the transition occurs for w ^ 1. The corresponding values for k\ 
in GeV^ for A = 0.3 and XI = 3 are also indicated. 



For small m- values the sum over k gives approximately exp (m/A), and thus 

a ~ exp(aoA/) ■ exp(M/A) = with A = aoA. (35) 

The transition between the regimes occurs now for a fixed ratio between u/l ^ 
A^/ao, which is of order 0.1 if A is around 0.3. The result is illustrated in fig. 
which shows how the expression in eq. (|33D interpolates between the DGLAP and 
BFKL limits in eqs. (|3^) and (^5|) for large and small /c_l.0 

The simple models in eqs. and (PB|), for fixed and running couplings respec- 
tively, interpolate smoothly between large and small k±-valnes. They contain the 
essential features of the CCFM and LDC models, and can therefore give an intu- 
itive picture of the transition between these two regimes. The transition occurs for 
a fixed ratio between log k\ and log l/x for fixed coupling, and between log(log k\) 
and logl/x for a running coupling. A more extensive discussion, including effects 
of non-leading terms in log 1/x and the properties of Laplace transforms, is found 



in ref. \ . 



^ In the framework of = 4 SUSY, where the corresponding coupHng is not running, the 
relations between DGLAP and BFKL equations has been found in [ |8^, |4| ^ for the first 
two orders of perturbation theory. These results can be useful for possible future study of the 
corresponding approximate relations in QCD. 
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4.4 The CCFM equation in the transverse coordinate rep- 
resentation 



It might be useful to consider the CCFM equation using the transverse coordinate 
(or impact parameter) [ ^ b conjugate to the transverse momentum k_i_ and to 
introduce the b dependent gluon distribution A{x, b, (f) defined by the Fourier- 
Bessel transform of A{x, k'j_, (f): 

_ POO 

A{x,b,q^)= k±dk±Jo{bk^)A{x,kl,q'^). (36) 
Jo 

where Jo{u) is the Bessel function. The advantage of this representation becomes 
particularly evident in the so called 'single -loop' approximation that corresponds to 
the replacement q/z ^ q and Aat^ — *■ 1. In this approximation the CCFM equation 
reduces to the conventional DGLAP evolution and eq.(p3D is diagonalized by the 
Fourier - Bessel transformation provided that we set the same argument /i^ ~ 
in both terms in this equation. The evolution equation for A{x, b, (f) reads: 

where 

P{z) = -^ 2 + z{l-z) + - (38) 

1 — z z 

In equation ( |37D the argument /i^ of the QCD coupling was set /i^ = if and 

A,(g^ Ql)=expi- ^ dzas{q')zP{z) (39) 

\ JQI q^ Jo J 

In the splitting function P{z) we have included the terms which are non-singular 
at 2; = or 2; = 1 in order to get a complete DGLAP evolution equation for 
the integrated gluon distribution. For the same reason we have also extended the 
definition of the Sudakov form factor by including the complete splitting function 
P{z) and not only its part which is singular at z = 1 (cf. eq. (p9|)). 
Equation ( pTl) can be solved exactly after going to the moment space and intro- 
ducing the moment function Auj{b,f) defined by: 

A^{b, f) = C dxx^A{x, b, f) (40) 
Jo 

The solution of the evolution equation for the moment function Auj{b, f) reads: 



X(&,r)=^"(fe)-exp 



/ ^asiq') dzzP{z){z^-'Ubq{l-z))-l) 
JQI q^ Jo 



(4 

where A^{b) is the Fourier- Bessel transform of the (input) unintegrated distribution 
at g = Qo- It should be noted that at 6 = 0, A^{b, (f) is proportional to the moment 
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function of the integrated gluon distribution g{x,q). The solution ( ^T] ) at 6 = 
reduces to the solution of the LO DGLAP equation for the moment function of 
the integrated gluon distribution. It should also be noted that the integrand in 
eq.(^lD is free from singularities at 2 = 1 and so we can set Qo = in the upper 
integration limit over dz. The formalism presented above is similar to that used 
for the description of the transverse momentum distributions in (for instance) the 



Drell-Yan process (see e.g. [ |8^, 

In order to obtain more insight into the structure of the unintegrated distribution 
which follows from the CCFM equation in the single loop approximation it is useful 
to adopt in eq. (|36|) the following approximation of the Bessel function: 

Join) - 6(1 - u) (42) 

which gives: 

Using the same approximation for the Bessel function in eq. (|4l|) we get, after some 
rearrangements, the following relation between unintegrated and integrated gluon 
distributions: 

A{x,kl,f) ^ ^[g(x,kl)A,{f,kl)] (44) 



Similar relation has also been used in refs. [ 0, |93| (see also eq.(p4D in the next 
section). 

In the general 'all loop' case, the non-Sudakov form-factor generates contributions 
which are no longer diagonal in the b space and so the merit of using this represen- 
tation is less apparent. However in the leading log(l/a;) approximation at small x 
the CCFM equation reduces to the BFKL equation with no scale dependence and 
the kernel of the BFKL equation in b space is the same as that in the (transverse) 
momentum space. The work to explore the CCFM equation in b space beyond the 
'single loop' and BFKL approximations is in progress. 



4.5 Available parameterizations of unintegrated gluon dis- 
tributions 

Depending on the approximations used, different unintegrated gluon distributions 
can be obtained (see eq.(|T5D). Here we describe some of the so far published 
ones and make some comparisons.^ We use the following classification scheme: 
xQ{x,k\) describes DGLAP type unintegrated gluon distributions, xJ-'{x,kj_) is 
used for pure BFKL and xA{x, k'j_, (f) stands for a CCFM type or any other type 
having two scales involved. 
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Parameterizations of the unintegrated gluon distributions described here are available as a 



FORTRAN code from tittp : //www . thep . lu . se/Smallx 
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Derivative of the integrated gluon density 

Ignoring the fact that the unintegrated density may depend on two scales we 
can invert eg . ([151) which gives the unintegrated gluon density G{x, k"^): 



(45) 



Here xg{x, fi'^) can be any of the parameterizations of the gluon density avail- 
able. 

GBW (Golec-Biernat Wiisthoff [0) 

Within the color-dipole approach of a simple parameterization of the 

unintegrated gluon density was obtained which successfully describes both 
inclusive and also diffractive scattering. The unintegrated gluon density is 
given by [ p^: 

H^,kl) = ^Rlix)klexp {-Rl{x)kl) (46) 
with (To = 29.12 mb, ag = 0.2 and 

1 /x\V2 



and Qo = 1 GeV, A = 0.277 and xq = 0.41 ■ 10"^, where the parameters 



correspond to the parameterization which includes charm [ R5 



RS (Ryskin, Shabelski [ |9f 
A parameterization of the unintegrated gluon density satisfying the BFKL 
equation but including also the non-singular parts of the gluon splitting func- 
tion as given by DGLAP is presented in [ 96]. The integrated gluon density 
xgi^Xjfi"^) is defined by: 

xgix, /i') = (j){x, kl)dkl (47) 

The function 0(x, /c^) is obtained by a solution of the evolution equation (in- 
cluding LO BFKL and the remaining part of the DGLAP splitting function, 
but without applying the kinematic constraint) in the perturbative region 
kj_ > Ql. The non-perturbative part of the gluon density is identified as the 
coUinear gluon density xg{x, Qq) at a small scale Qq ~ 4 GeV^: 



:F{x,kl) = { ^ (45 



An explicit form of the parameterization of the function together with the 



fitted parameters are given in [ 96 
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KMS (Kwiecinski, Martin, Stasto 



In the approach of KMS [ ^ three modifications to the BFKL equation 
are introduced in order to extend its vahdity to cover the full range in x. 
Firstly, a term describing the leading order DGLAP evolution is added. The 
inclusion of this term has the effect of softening the small x rise and to 
change the overall normalization. Secondly, the kinematic constraint eq.(|5^ 
is applied, which has its origin in the requirement that the virtuality of the 
exchanged gluon is dominated by its transverse momentum. Thirdly, the 
BFKL equation is solved only in the region of k± > k_i_o, whereas the non- 
perturbative region is provided by the collinear gluon density 
the scale k±o. The border between the perturbative and non-perturbative 
regions has been taken to be k±o = 1 GeV^. 

Also a term which allows the quarks to contribute to the evolution of the 
gluon has been introduced. The single quark contribution is controlled 
through the g ^ qq splitting. The input gluon and quark distributions 
to this unified DGLAP-BFKL evolution equation are determined by a fit 
to HERA F2 data. The unintegrated gluon density J-'{x, k'j_) is still only a 
function of x and k'j_, which strictly satisfies: 

X9{x, Q') = J -^f{x, kl) = j dkl:F{x, kl). (49) 
JB (J. Bliimlein [|8[) 

In the approach of JB [ ^ the general fc^-factorization formula: 
(t(x, yu^) = j dk']_a{x, /c^, n^) ® A{x, k]_, fx^) 



is rewritten in the form \ 36 



a{x, yu^) = o-°(a;, n^)^g{x, J dk^ (^a{x, k\, n"^) - a°{x, n^)'^(S)A{x, k\, fi^), 

(50) 

which gives the relation: 

g{x, fi'^) = / dk']_A{x, k']_, ^'^). 
Jo 



In the case of eq.(|SD|) the k± dependent gluon distribution satisfying the 
BFKL equation can be represented as the convolution of the integrated 
gluon density xg{x,^'^) (for example GRV [ 15^) and a universal function 

B{x,kl,^l')■. 

A{x, kl, /i^) = B{z, kl, fi') - g{-, /i^) dz, (51) 

Jx Z Z 

The universal function B{x, kl, /i^) can be represented as a series (see [ p^): 

N=oo 

J2 dili-i{...), 

i=l 
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where = Jj if > fi"^ and Jj = Jj if k"^ < n^, respectively and Jj and Jj are 
Bessel functions for real and imaginary arguments. The series comes from 
the expansion of the BFKL anomalous dimension with respect to as- The 
first term of the above expansion explicitely describe BFKL dynamics in the 
double-logarithmic approximation: 



B{z,kl,ix' 



B{z,kl,fi') 



as 
zkl 



Jo(2v/aslog(l/;^)log(/iVA;i)), 



/o(2Jaslog(lA)log(fci//x2)), 



kl < /^^ 



kl > fi^, 



(52) 



(53) 



where Jo and Iq are the standard Bessel functions (for real and imaginary 
arguments, respectively) and as = Sa^/Tr is connected to the pomeron inter- 
cept A of the BFKL equation in LO cull = Q;s41og2. An expression for cjnll 
in NLO is given in [ ^ ^ : cjnll = tts4 log 2 — Na^. Since the equations be- 
have infrared finite no singularities will appear in the collinear approximation 
for small k i . 



KMR (Kimber, Martin, Ryskin [ [TOOi ) 

In KMR [ |100|| the dependence of the unintegrated gluon distribution on the 
two scales k± and fi was investigated: the scale fi plays a dual role, it acts as 
the factorization scale and also controls the angular ordering of the partons 
emitted in the evolution. Already in the case of DGLAP evolution, uninte- 
grated parton distributions as a function of the two scales were obtained by 
investigating separately the real and virtual contributions to the evolution. 
In DGLAP the unintegrated gluon density A{x, k\, fi^) can be written as: 



xA{x,k']_, 1^) = T(/c^,/i^ 



1 asik^) /-(i-^) 



T{kl,fi' 



kl 271 

dxg{x, /i^ 



P{z)x'g{x,kx)dz (54) 



where Tiki, ji^ 



is the Sudakov form factor, resumming the virtual correc- 



tions: 



T{k 



exp 



2 2tt dk? 



P{z')dz' 



The structure of eq.(^4D is similar to the differential form of the CCFM equa- 
tion in eq. (p3D , but with essential differences in the Sudakov form factor as 
well as in the scale argument in eq . (|5^) . The splitting function P{z) in eq. (|5^ 
is now taken from the single scale evolution of the unified DGLAP-BFKL ex- 
pression discussed before in KMS [ . As in KMS the unintegrated gluon 

1 GeV^. An extrapolation to 
such, that: 



density /(x, k^, ^'^) starts only for k\ > k 



±0 



cover the whole range in kj_ has been performed [ 101 



— P 



if A; J 
if A;, 



< k±o 
> k±o 



(55) 
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with xg{x,k'j_Q) being the integrated MRST [ [I02|| gluon density function. 
The unintegrated gluon density xA{x, k'^, fj,'^) therefore is normahzed to the 
MRST function when integrated up to kj_Q. 



JS (Jung, Salam [[T^, |54 



The CCFM evolution equations have been solved numerically [ |1J, ^4 



us- 



ing a Monte Carlo method. The complete partonic evolution was generated 



according to eqs. (p!q-|27|) folded with the off-shell matrix elements for boson 
gluon fusion. According to the CCFM evolution equation, the emission of 
partons during the initial cascade is only allowed in an angular-ordered region 
of phase space. The maximum allowed angle H for any gluon emission sets 
the scale q for the evolution and is defined by the hard scattering quark box, 
which connects the exchanged gluon to the virtual photon (see section [4.2|) . 

The free parameters of the starting gluon distribution were fitted to the 
structure function F2{x,Q'^) in the range x < 10~^ and > 5 GeV^ as 
described in [ resulting in the CCFM unintegrated gluon distribution 
xA{x, kj_, q^). 

In the following we show comparisons between the different unintegrated gluon dis- 
tributions. The JS unintegrated gluon density serves as our benchmark, because 
calculations based on this have shown good agreement with experimental mea- 
surements, both a HERA and at the TEVATRON. In Fig. |^ we show the gluon 
distributions at /i = g = 10 GeV obtained from the different BFKL sets [KMR, JB) 
as a function of x for different values of k'j_ and as a function of kj_ for different val- 
ues of X together with the JS unintegrated gluon density. It is interesting to note, 
that the JS gluon is less steep at small x compared to the others. However, after 
convolution of the JS" gluon with the off-shell matrix element, the scaling violations 
of F2{x,Q'^) and the rise of F2 towards small x is reproduced, as shown in [ [l^. 



Fig. 4 therein]. In the lower part of Fig. |^, the effect of the evolution scale on the 
distribution in kj_ is shown: The JS and KMR sets both consider angular ordering 
for the evolution, whereas in JB the evolution in k_i_ is decoupled from the evolu- 
tion in fi, resulting in a larger kj_ tail. JS also includes exact energy-momentum 
conservation in each splitting which further suppresses the high-fc^ tail. Fig. ^ also 
shows, that the small k"^ region, which essentially drives the total cross sections, 
is very different in different parameterizations. 

In Fig. 1^ we show a comparison of the gluon density distribution obtained from the 
derivative method (using GRV98 LO) and KMS and compare it to the JS gluon 
at g = 10 GeV. The KMS and "derivative of GRV" give very similar unintegrated 
gluon distributions, which is a result of the strict relation between the collinear 
(integrated) and unintegrated gluon distributions, that was used in KMS. One 
also has to note, that KMS and GRV98 use very similar data sets from HERA 
for their fits. In Fig. ^ we show a comparison of the unintegrated gluon density 
parameterizations from GBW and RS and compare it to the JS gluon at g = 10 
GeV. The RS set is shown only for historical reasons, since it was one of the first 
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Figure 7: The unintegrated, two scale dependent gluon distributions at q = jj, = 
10 GeV as a function of x for different values of k'j_ (upper) and as a function of 
k'j_ for different values of x (lower): JS [^^, (solid line), KMR [ \10(\j (dashed 
line) and JB (dotted line) 



unintegrated gluon distributions available. The (75 unintegrated gluon density, 
although successful in describing inclusive and diffractive total cross sections at 
HERA, is suppressed for large k± values (see Fig. which is understandable, since 
large k± values can only originate from parton evolution, but were not treated in 
GBW. However, it is interesting to compare it also to more exclusive measurements. 
One can also see, that the Gi^VK gluon decreases for /c^ — s> 0. 
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Figure 8: The unintegrated, one scale gluon distributions as a function of x for dif- 
ferent values ofk\ (upper) and as a function of for different values of x (lower): 
KMS /[g^/ (dashed line) and derivative GRV LO (dotted line) here compared to the 



two scale gluon distribution of JS [^^, [7^/ (solid line) at q = jj, = W GeV (same 
as in Fig.^ 



In Table ^ we present cross sections calculated with four of the different uninte- 
grated gluon distributions discussed above for heavy flavor production and inclusive 
deep inelastic scattering at HERA energies {^/s = 300 GeV) as well as the total 
bottom cross section at the TEVATRON. The benchmark is again the JS param- 
eterization, which is able to describe the corresponding measurements at HERA 
and at the TEVATRON reasonably well. Large variations in the predicted cross 



28 




Figure 9: The unintegrated, one scale gluon distributions as a function of x for 
different values of k\ (upper) and as a function of k'j_ for different values of x 
(lower): GBW (dashed line) and RS f^B{] (dotted line) here compared to the 
two scale gluon distribution of JS 



J^/ (solid line) at q = jj, = 10 GeV (same 



as m 



FtgW 



sections are observed. It is clear that the parameterizations of the unintegrated 
gluon distributions are very poorly constrained, both theoretically and experi- 
mentally. The differences in definition makes it difficult to go beyond a purely 
qualitative comparison between the different parameterizations. Also, since even 
the integrated gluon density is only indirectly constrained by a fit to F2{x, Q^), it 
may be necessary to look at less inclusive quantities to get a good handle on the 
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(T[nb] 




JS 


KMR 


JB 


GBW 


ep e'ccX (Q2 < i GeV^) 


696.8 


412.8 


741.3 


735.3 


ep e'ccX (Q2 > i GeV^) 


80.2 


47.4 


87.7 


82.0 


ep e'bbX (g2 < 1 GeV2) 


5.36 


2.78 


4.38 


5.64 


ep e'X (Q2 > 1 GeV2) 


838.7 


610.0 


6550.0 


4134.0 


ep e'X (Q2 > 5 GeV^) 


212.8 


127.2 


585.1 


564.0 


pp bbX (v/i = 1800 GeV) 


88100. 


27489. 


78934. 


65990. 



Table 3: Cross sections for different processes at HERA and TEVATRON using dif- 
ferent parameterizations of the unintegrated gluon distributions. The cross sections 
are calculated with the Cascade /|53|/ Monte Carlo generator. In all cases the one- 
loop as{p^) is used with p? = p\-\-m'^, where p± is the transverse momentum of one 
of the quarks in the partonic center- of-mass frame and nig = 0.140, 1.5,4.75 CeV 
is mass of the light, charm and bottom quarks. 



unintegrated gluon. 



4.6 Beyond leading logarithms 

In this section the attempts to go beyond leading order are described and summa- 
rized. First, general aspects of next-to-leading effects are discussed and then the 
CCFM approach is critically considered. 

We write the momenta of the t-channel gluons (in Sudakov representation) as 
ki = xfPi -\- x~ P2 + k^i and the emitted gluons as Pi = Vi{Pi + ^4^2) + P±i 
(eq.(^)). The discussion in this section will be based on the strong ordering limit, 
that is x~^, angles and x~ are strongly ordered, which means that factors of (1 — z) 
can be safely ignored since they are of 0{1) for z ^ 0. In the high energy limit 
with strong ordering, one can talk of x~^ ordering of exchanged gluons or of angular 
ordering of the emitted gluons. The variable C, (related to the angle of the emitted 
gluons) can be written (in this approximation) as: 



6 



X. 



(56) 



i-l 



with Vi ~ x'l_i for 2 ^ (eq.([T^)). The last expression in the above equation is 
obtained from: ^ 



X, 



i-l 
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4.6.1 General next-to-leading order investigations 



It is known from many contexts of QCD that for reasonably accurate predictions it 
is necessary (at the very least) to go to next-to-leading order. In what follows, one 
has to remember that the role of next-to-leading corrections in the fc^-factorization 
approach is very different to the ones in the collinear approach, since part of the 
standard NL corrections are already included at LO level in /c_L-factorization, as 
was discussed in chapter ||. 

For a next-to-leading logarithmic (NLL) calculation of a cross section at high en- 
ergies, there are two main ingredients. One is the NLL corrections to the 'kernel' 
of the BFKL equation, generating terms as{aslogl/x)"'. This part should be in- 
dependent of the process under consideration. The second part is the correction 
to the impact factors (off-shell matrix elements) at either end, and is the source of 
the process dependence in the NLL corrections. Processes of interest include 7*7* 
scattering in electron positron annihilation, forward jets at HERA and Mueller- 
Navelet jets at hadron-hadron colliders. In order to describe these processes at 
NLO level of accuracy, one needs the photon impact factor and the jet vertex. For 
both elements NLO calculations are on the way: for the photon impact factor the 
main are the off-shell matrix elements described in section |3]^ (refs. [ ^ [70[]), 
and the quark induced jet vertex has recently been computed in [ |TU^. 

One of the major developments in past years has been the completion of the calcu- 
lation of the NLL corrections to the BFKL kernel. This was a significant enterprise, 
lasting almost a decade [ Once the various contributions have been as- 



sembled, the final result can be given in a fairly compact form [ |3^, ^ . It can 
be summarized through the following relation between the next-to-leading BFKL 
power, cunll, and the leading power cull = as41og2: 

c^NLL = «s4 log u;^^{l - 6.2a,) . (58) 

Substituting a typical value for a,, say 0.2, one finds a negative power — so rather 
than improving the accuracy of the predictions, the NLL corrections seem to lead to 
nonsensical answers. A more detailed analysis suggests that for problems involving 
two substantially different transverse scales, the inclusion of the NLL corrections 



leads to an even worse problem, namely negative cross sections [ 104 |. So at first 
sight it seems that the perturbation series for BFKL physics is simply too poorly 
convergent for it to be of any practical use. 

Despite this, there are indications that ways exist of using the NLL corrections 
for phenomenological purposes (7*7*, forward jets or TEVATRON jets with large 
rapidity separation might be examples for this). This is because it is possible to 
identify a well-defined physical origin for large parts of the NLL terms. These 
parts can then be calculated at all orders, and the remaining pieces are then a 
much smaller NLL correction. 

A clue as to the origin of the large corrections can be obtained by examining their 
structure in the collinear limit (where one transverse scale is much larger than 
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the other). From DGLAP we think we understand the origin of all terms involving 
{as log Q^)" - they are associated with strong orderings in k±. So e.g. in the BFKL 
NLO corrections (terms oc a^) we have a piece with (a^ log^ Q^) which is something 
we already know about from DGLAP. Terms with this "collinear" enhancement (a 
number of powers of log Q) turn out to be responsible for over 90% of the NLO 
corrections to BFKL, and so are the reason for the large size of these corrections. 
But since their origin is just DGLAP physics, which we know well, we can also 
predict the terms that arise at NNLO (a^ bg^ Q^) etc. and resum them. 

Suppose one wishes to calculate the high-energy behavior of a Green's function 
with a squared center of mass energy s and transverse scales Q,Qo at the two ends 
of the chain. It is convenient to examine this in Mellin transform space, with u 
conjugate to s/{QQq) and 7 conjugate to the squared transverse momentum ratio 
Q'^/Qq- One can then write the BFKL kernel as 

ctsNc 

= as ixo + asXi H ) , «s = • (59) 

TT 

For small 7 (corresponding to a large ratio of transverse momenta) the leading part 
of the BFKL kernel goes as 

Xo ^ - + 0(f) • (60) 

7 

In the same region the NLL corrections behave as 

+ (61) 

The extra divergences at small 7 can be quite easily understood because each power 
of 1/7 (after inverse Mellin transform) corresponds to a logarithm of transverse 
momentum (logQ). So for example the term with 1/7^, given that it multiplies 
Og, corresponds to a contribution (a^ log Q^/ Q^)"^, and so looks like a standard term 
from DGLAP evolution. One can verify the coefficient that would be expected from 
leading-log DGLAP evolution and it comes out as exactly —11/12. 

The term proportional to 1/7^ is slightly more subtle because it looks super-leading 
compared to DGLAP. Recalling that it is multiplied by ag, we can see that it is 
related (after a Mellin transform) to a term which we can write as (a^ log Q^/ Ql) x 
{as log^ Q^/Qo)' i-^- ^ normal DGLAP term, multiplied by a double log of Q^. Such 
double logs of will be discussed in more detail in the next subsection. Essentially, 
their presence or absence depends on whether we resum high energy logs oix~^, x~ 
or -y/x+Tx" (or equivalently Q^/s, Qq/s or QQo/s). The small-7 part of the NLO 
BFKL kernel has been shown above with the high-energy logs defined in terms of 



/x — converting to the natural evolution variable for DGLAP evolution, 
the 1/7^ pole disappears, and one is left with just the 1/7^ term. In other words 
the small-7 behavior is entirely constrained by known results from DGLAP. 

It is important to note though that the 1/7^ term is not purely an artefact in 
that there is an analogous triple pole around 7 = 1, associated with transverse 
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logarithms for Qo ^ Q, and there is no convention for the definition of the high- 
energy logs which allows one to get rid of the triple poles simultaneously at 7 = 
and 7 = 1. 

Given that these coUinearly enhanced contributions have a simple physical origin, 
namely the strong ordering in transverse momenta, they can be quite straightfor- 
wardly calculated not just to NLL, but to all orders. For example the generalization 
of the 1/7^ and 1/(1 — 7)^ terms can be obtained by implementing the kinematic 
constraint discussed above. The generalization of the 1/7^ and 1/(1 — 7)^ terms 
comes from the leading-log DGLAP kernel and from the running of the coupling. 

Of course an understanding of the structure of the kernel around 7 = does not 
formally tell us anything about 7 = 1/2, which is the region of interest for small- 
X predictions. What is remarkable though is that a collinear approximation, i.e. 
taking just the known poles at 7 = and 7 = 1, 

coll _ 1 1 11 11 

" ~ 2(1 - 7)3 " 12y ~ 6(1 - 7)2 ' ^^^^ 

approximates the full next-to-leading corrections to better than 93% over the whole 
range < 7 < 1. The above equation is just the expression in '7'-space (the Mellin 
transform of Q space) of terms that we know from DGLAP. 

This suggests that a practical approach to 'improving' NLL BFKL might just be to 
include the higher-order collinear terms at all orders. This has been done in [ p.U5| , 
106| , |107|| , and is found to considerably improve the stability of the predictions for u, 



giving a phenomenologically sensible value of cu in the range 0.25 to 0.3 for ~ 0.2. 
These kinds of improvements are necessary whenever exact NLO corrections are 
used, because without them the result is very unstable (large renormalization scale 
uncertainty, negative power u for ag > 0.15). The same approach can be extended 
to the calculation of Green's functions and physical cross sections, and work in this 
direction is currently in progress [ |108|| . 

There have been other approaches to supplementing the NLL corrections with 
higher-order contributions, notably performing a BLM (Brodsky Lepage Macken- 
zie) change of the scale for [ |109[ | and enforcing rapidity separations between 
emissions [ |110|| . The basic idea behind the BLM scale choice is that a large part 
of the NLO corrections come about because in the LO calculation the scale choice 
for as was unreasonable. This is reflected in the NLO corrections by a large piece 
proportional to /3o, which is one of the parameters of the /3 function controlling 
the renormalization scale dependence of a^. However, in a NLO calculation one 
obtains coefficients of Uf, and of C^, but not of /3q as such, so one can't identify 
what the actual coefficient of Pq is. BLM identifies the "amount" of jSo with the 
entire nj part, effectively guessing from the coefficient of n/ the coefficient of (3o, 
and then saying that it should all be re-absorbed into a different choice of scale 
for ttj, in the LO piece. When the authors of [ |109|| tried this for NLO BFKL they 



found that it made the calculations marginally worse (even more unstable), unless 
they chose a very specific renormalization scheme. They argue that this scheme. 
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which appears in certain other contexts involving three-gluon couphngs (though 
not particularly in a small-x limit) might be a more natural choice when applying 
the BLM procedure, and find a value for uj of around 0.15. 

The approach which enforces rapidity separations is motivated by the observation 
that the approximations in the LL calculation for BFKL are satisfied only when 
there is a large gap in rapidity between successive emissions, but that one never- 
theless integrates over all possible separations between emissions (up to the limit 
where emissions have the same rapidity), which is plainly in violation of the initial 
approximations. One can choose instead to integrate up to some (a priori arbi- 
trary) minimum rapidity separation, Ar/. This modifies one's prediction at higher 
orders, and in particular can partially mimic the NLL corrections. However its 
explicit merging with the NLL corrections by Schmidt [ |110|| revealed a significant 
remaining instability with respect to variations of Ar/. There is evidence though 
that when combined with the coUinear resummation described above, much of 



the instability with respect to the variation of A?7 disappears [ Essentially 



both methods (coUinear enhancement, rapidity separation) cut out certain regions 
of phase space. With the "coUinear" approach the piece cut out is well defined, 
whereas with the rapidity cut, the amount cut out depends on Ar/. But when the 
two methods are put together, a large fraction of the phase space that would be 
removed by a rapidity cut has already been removed by the coUinearly enhanced 
terms, so the effect of rapidity cut is much reduced. 



4.6.2 Anomalous dimensions. 

So far the discussion has been most relevant to processes with two similar hard 
scales. In the case of re-summed small- a; splitting functions (anomalous dimen- 
sions) in addition to the issues discussed above, there is a further subtlety related 
to iteration of the running of the coupling. It has been discussed with comple- 
mentary methods but similar conclusions in [ |1(J7| , |11(J , TVl], [ ITS, 114, 115| and 



|116|| . Essentially the result is that for a small-z splitting function at scale Q^, the 
typical scale in the evolution is considerably higher than 0. Similar results have 
been found before also in [ |12g, [121], Q g in the framework of DGLAP, and 



in [ p.09|| in the framework of BFKL. The splitting function is only independent of 
at leading order. At NLO it is no longer scale invariant because the NLO piece 
comes in with a different relative weight according to ^^(Q^), as do NNLO etc. 
When deducing the small-x enhanced part of the splitting function then the NLO, 
NNLO, NNNLO (in the DGLAP hierarchy, <+'"log"g2 ^ N'^LO) pieces which 
are small-x enhanced need to be determined at all orders. And the final 'splitting 
function' then depends on because the relative weights of all these pieces is 
different according to the value. The inclusion of all these terms causes the 



The effective fiigh scale could be responsible (see discussions in [ 117, 118, 119|) for the good 
agreement between the experimental data for the F2 structure function and the pcrturbative 
estimations in the small x range. 
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Pgg{z) splitting function to grow as a power l/z^^'^" at small z. Initial calculations 
of this dependence led to a large value for uic and consequent incompatibility with 
the data (see e.g. [ |124|| ). However these predictions are modified by two classes 
of large effects: NLL contributions, similar to those discussed above, and also cor- 
rections associated with the fact that the effective scale of relevant to small-2; 
splitting functions is substantially larger than Q^. After an involved calculation 
one finds that, parametrically, this implies corrections to ujc of order a^^^, aj^'^, . . . 
(which need to be resummed) and leads to a significant reduction in the small-x 
power of the splitting function compared to the fixed-coupling CclSG, clS well as in- 
creased stability with respect to the NLL corrections. This means that the power 
ujc for the small- 2; dependence of the splitting function, l/z^^'^", is smaller than 
one would expect (i.e. it's not the same Lj{as{Q'^)) that is calculated for a process 
such as 7*(Q)7*(Q)) and furthermore it sets in only at very small z (below 10^^). 

An approach which also involves some of the elements discussed above has been 



proposed and used in [ |125| , |126|| to study scaling violations in the HERA F2 data 



|127|| . In their approach the authors make use of a coUinear resummation around 
7 = 0, however argue that this cannot also be applied to 7 = 1, which implies that 
it is not possible to predict the height of the minimum of the characteristic function. 
Accordingly the asymptotic power A of the small-a; splitting function ~ x"'*' must 
be fitted to the data. Within this approach the authors find an improved agreement 
with the F2 data compared to a pure NLO DGLAP fit. The optimal value for A 
depends on the details of how the resummation is implemented, in one case being 
negative (A ~ —0.25), while in the other it is of the order of A ^ 0.2 which would 
also be consistent with other resummation approaches. 



4.6.3 CCFM and problems beyond leading logs 

The CCFM evolution involves three partonic variables, the light-cone momen- 
tum fraction, the emission angle and the transverse momentum. Looking at the 
case where a forward jet in DIS has a p\ much smaller than Q^, the natural 
variable when going up in scale from the proton towards the virtual photon is 
the positive light-cone momentum fraction, = Q'^/{sy) ~ Xbj, giving rise to 

the standard double-logs, (^as log ^ log , coming from integrals of the form 

1^2 dq^ / dz / z. Naively going from the photon side with the negative light- 
cone momentum fraction (compare eq. (|5^) ) x~ = p]_/{sy) = x^p]_/Q'^ will give us 
so-called illegal double-logs of Q"^ /p\ 

as log — log — = as log — log — + as log — (63) 
V P± ^~ ) \ Pi a;+ Pi) 

The wording illegal comes from the observation, that the renormalization group 
equations (DGLAP) tell us what classes of logs we can expect (they are first order 
differential equations so we expect only single logs of Q'^). When we get a double 
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log of then this is illegal because it cannot be compatible with the renormaliza- 
tion group. Subtleties, of course, arise because different ways of writing the same 
expression can lead to the double log being visible or not, and that is what the 
above eq.(^) describes. CCFM uses the angle (or more specifically ^) as evolution 
variable. Using = then the ordinary coUinear evolution with increasing 

p± and automatically decreasing results in angular ordering or ordering in ^/^. 
If on the other hand the p±^s decrease then a decrease in does not automatically 
imply an increase in angle - so angular ordering becomes a stronger condition. It 
is the angular ordering ~ ^Jx~^/x^ which is responsible for the illegal double 
logs. 

In the original formulation of CCFM a kinematic constraint was noted: 



which was derived already in [ 128 , eq.(2.11)] and in [ 129 | in a frame where 
parton i is along the z axis. If the transverse momenta k±i are to be the dominant 
contributions to the virtualities /cf (as assumed in the derivation of the CCFM 
equation) then the so-called kinematical or consistency constraint is obtained: 

/cL > z+,kl,_,. (64) 

To merely have a self-consistent evolution equation (no infrared unsafety, proper 
real- virtual cancellations), which gives correct cross sections and final-state distri- 
butions at leading logarithmic accuracy, the consistency constraint is not needed 
in the CCFM equation. The way the consistency constraint was originally writ- 
ten into the CCFM equation should not be used because it messes up real-virtual 



cancellations as discussed below in section 4.6.4 



4.6.4 NLL corrections and event generators 

There are also aspects in which the NLL corrections may be able to help guide 
the construction of event generators. They give a hint about the kind of physical 
effects which are likely to be important (splitting functions, kinematic constraint). 
There are also situations when building an event generator in which one is faced 
with two possible strategies, and sometimes the NLL corrections give clues as to 
the direction to be followed. 

One example relates to the running of the coupling. A priori it is not necessarily 
clear what scale to use and for simplicity, in various phenomenological contexts, 
the scale has often been chosen to be k^- (as related to the branching (i — 1) ^ i 
in figure 0). However, in the full NLO calculation for BFKL a term appears 

-ai/5oln%. (65) 

The interpretation is very simple: the right scale in the LO kernel is simply the 
emitted p± (because if included in the LO kernel and then expanded in powers of 
a<j, one obtains precisely such a term at NLO). 
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But knowledge of the 'ingredients' needed for an event generator is not always suf- 
ficient. One generally has a basic equation for the branching, such as the CCFM 
equation, but inserting the relevant corrections such as to obtain the correct struc- 
ture of logs at the end is often not trivial, because the full set of logs comes out of 
the iteration of the branching. Furthermore some care is often needed to maintain 
proper cancellation of real and virtual corrections. 

For example in normal BFKL the correct cubic poles around 7 = 1 can be obtained 
by inserting the requirement kf > Zikf_i for each branching. This corresponds to 
symmetrizing evolution up in transverse scale (with as the evolution variable) 
and evolution down in transverse scale (with x~). 

However implementing this constraint directly into CCFM leads to the wrong coef- 
ficient for the 1/(1 — 7)'^ term — in other words the symmetrization is not properly 
accomplished. One can partially solve the problem by modifying the virtual cor- 
rections (non-Sudakov) form factor 

In A„3(., g^ kl) = -f^f' ^as(/)e(g'^ - z'V)eikl - g'^)e - z] , 

(66) 

where the last O function accounts for the kinematic constraint in the virtual 
corrections. This will not make the cross section exactly symmetric, but it will 
ensure that for evolution downwards in scale the evolution variable does at least 
correspond roughly to x~ . On the other hand as far the final state structure is 



concerned, this 'hack' [ [L3CI|| cuts out certain regions of phase space that should be 
left in. An a priori way of determining whether the 'hacks' are any good should 
be to examine how they compare to exact NLO calculations (but that means 
determining their expansion to 0{al), which isn't necessarily an easy task). 

Another source of next-to-leading-log corrections is the gluon splitting function. 
At very large energies, the 1/z term in Pgg, included in BFKL and CCFM, will 
certainly be dominant. However, the question is whether the treatment of just this 
term is sufficient at energies available at present colliders. 

The effect of the new small x is best seen in forward jet production in deep inelastic 
scattering, where the contribution from typical DGLAP dynamics is suppressed. 
In such a process, using a Monte Carlo simulation, the distribution of 2- values can 
be studied, and the validity of the small x approximation can be checked. 

The basic event selection criteria for forward jet production at HERA are given in 



Fig. 10 



The criterion E^/Q^ is essential to suppress the DGLAP contribution within a typ- 
ical range in x of 10~^ < x < 10~^. The evolution takes place from the large Xjet 
down to the small x with a typical range at HERA energies of Ax = x/xjet > 0.01. 
In order to justify the use of an evolution equation (instead of a fixed order cal- 
culation) one would require at least 2 or more gluon emissions during the evolu- 
tion. To roughly estimate the energy fractions Zi of 3 gluon emissions between 
10~^ < X < 10^^, one can assume that all gluons carry equal energies. Then the 
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g2 > 10 GeV2 
ETjet > 5 GeV 
Vjet < 2.6 
Xjet > 0.036 
0.5 < E^Q^ < 2 




10 



10 



Figure 10: The values of the splitting variable z for events satisfying the shown 
forward jet criteria, with 6 = 7° at HERA energies. 

range of Ax ~ 0.01 results in 2; ~ 0.2, which is far from being in the very small z 
region, where the BFKL or CCFM approximations (treating only the 1/z terms in 
the gluon splitting function) are expected to be appropriate. In Fig. ^ we show 
the values of the splitting variable z in events satisfying the forward jet criteria at 
HERA energies obtained from the Monte Carlo generator Cascade [[1^. Since 



the values of the splitting variable z are indeed extending into the medium and 
large z region (the majority has z > 0.1), it is questionable, whether the BFKL 
evolution equations, including only the 1/z part of the gluon splitting function, 
are applicable. CCFM also has a 1/(1 — z) term, which means that in any case 
one expects quite a few emissions for 2; 1, but still, the medium z range is 
questionable. 

The implementation of the full DGLAP splitting function into CCFM is problem- 
atic. Naively one would simply replace 

- 2 + 2^(1 - z) , (67) 

1 — z 1 — z 

in the CCFM splitting function. However doing this leads to negative branching 
probabilities, so one needs to be a little more 'subtle'. 

For example a positive definite branching probability can be obtained by making 
the following replacements 

- - ^ + Bz{l - z) , J_ ^ ^ + (1 - B)z{l - z) , (68) 
z z 1 — z 1 — z 
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resulting in: 

P{z,q,k) = as{^kl)(^^^^^ + {l-B)z{l-z)jAns{z,q,k) (69) 

+as ((1 - zfq') + Bzil - z)^ 

where 5 is a parameter to be chosen arbitrarily between and I.Asa consequence 
of the replacement, the Sudakov form factor will change to: 

log A,(g^gg) ^_ IJ_ ^^^^ (^2^1 _ ^)2^) ^ + (1 _ 5)^(1 _ ^) 

(70) 

and also the non-Sudakov form factor needs to be replaced by: 

log A.. = -as {kl) £dz' + (1 - B)zil - z)) J ^Q{k - q'Mq' - z'q) 

^ (71) 

which can be rewritten 

log = -a, {kl) dz' {^-^ + (1 - B)z{l - z)^ 2 In ^Q{k - z'q) (72) 

The proposed changes to the original CCFM splitting function are 'hacks'. There 
are reasonable arguments, to change the scale in to q^{\ — zY everywhere, 
to include the full splitting function instead of only the singular parts and also 
to impose the kinematic constraint, with its consequences. These 'hacks' are of 
course constrained by the requirement that the NLL corrections corresponding to 
the given branching equation are similar to the true NLL corrections. However, 
this will not be the real final solution and perhaps it is a need for a better starting 
point, embodying both angular ordering and symmetry right from the start. 



5 Generators for small- cc evolution 

Three different Monte Carlo event generators are available, specifically devoted to 



small X processes (in chronological order): Smallx [ LDCMC [ ^ [51 



52, 53 1 and Cascade [14, 54, M 



5.1 Smallx 



Smallx [ 48, 49 is a Monte Carlo event program which generates events on parton 
level at small x in ep scattering. It uses the CCFM [41, 42 , 43 , 44 1 evolution 
equation for the initial state cascade convoluted with the /c^-factorized off-shell 
matrix elements of [ |3^ for light and heavy quark pair production. Modifications 
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of the original Smallx [ HBI, HH] version concerning the non-Sudakov form factor 



Ans were necessary to obtain a reasonable description of the structure function 
F2{x, Q"^) as well as hadronic final states like the forward jets at HERA [ |131| , |54 



In Smallx the initial state gluon cascade is generated in a forward evolution 
approach. The gluon evolution starts from the proton side with an initial gluon 
distribution according to (including a Gaussian intrinsic distribution around 
fco): 

XoG,{x,, kl,) =N-{1- xo)' ■ exp {-kljkl) (73) 

with being a normalization constant. A gluon with momentum fraction Xj and 
transverse momentum k±i is allowed to branch into a virtual (t-channel) gluon 
with momentum fcj+i and a final state gluon with momentum p^+i according to 



the CCFM splitting function [ ^ This procedure is repeated until the 



next emitted gluon would violate the angular bound q given by the matrix element. 

Smallx generates the full parton level structure, but due to the complicated struc- 
ture of the initial state branchings and the phase space of the matrix element, a 
weight is associated with each event. Although Smallx produces weighted events 
and therefore is inefficient for generating specific exclusive signatures, it can be 
used for the CCFM evolution and to calculate the inclusive structure function 
F2{x,Q'^). By performing a fit to measurements of the inclusive structure func- 
tion F2{x,Q'^), Smallx can be used to determine the unintegrated CCFM gluon 
density xA{x, k'j_, (f) in a grid in x, k^ and (f (with q being the scaled maximum 
angle allowed for any emission with (f = x^_i'Es). This numerical representation 
can be used for any other calculation. It is advantageous, that the full parton 
level is generated during the evolution, since this information can be used for com- 
parison with other event generators adopting a more efficient backward evolution 
approach [ |53, |55|. 



5.2 Cascade 



Cascade [ |T^, ^ is a full hadron level event generator which uses k±- 
factorization of the cross section into an off-shell matrix element and an uninte- 
grated gluon density function. The initial state cascade is generated in a backward 
evolution approach, which was necessary for an efficient event generation. The 
Lund string fragmentation package Jetset/Pythia [ |132|| is used for hadroniza- 



tion. Cascade can be used in ep, 'jp and also pp processes. 

The hard scattering process is calculated using the off-shell matrix elements given 
in [ ^ for light and heavy quark pair production, or ^yg ^ J/ipg convoluted 
with the unintegrated gluon density xA{x, k"^, (f) . It could be shown in [ |5^, |5^ 

5^ , |133| , |134|| is possible for small x 
ordering in the initial state cascade. 



that a backward evolution approach [ |I^ 
processes, which are not restricted to strong fcj 
This was the main ingredient for the development of a time efficient Monte Carlo 
event generator. 
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The backward evolution starts from the hard scattering process and evolves the 
partons backwards towards the proton. This approach is much more efficient, com- 
pared to the strategy of Smallx which uses a forward evolution approach. How- 
ever, it requires the unintegrated gluon density to be determined beforehand. The 
CCFM unintegrated gluon distribution has been determined from a Monte Carlo 
solution of the CCFM evolution equation which has been fitted to the measured 
structure function F2{x,Q'^). It is advantageous to have the unintegrated gluon 
density determined in a Monte Carlo approach with full control over the partonic 
state, as available in Smallx, since this allowed to show [ |14| that the backward 
evolution produces identical results to the forward evolution approach on the par- 
ton level. The proof of equivalence of the forward and backward evolution is a 
unique feature of CCFM. For DGLAP this has only been shown at at an inclusive 
level, while at the exclusive level, complications can arise related to differences in 
the treatment of angular ordering between forward and backward evolution. 

The program code is available from http://www.quark.lu.se/~hannes/cascade 



5.3 LDCMC 

Since the LDC model is inherently forward-backward symmetric, it is natural to 
design a Monte Carlo in the same way. In the LDCMC program all emissions are 
generated in one go from an incoming, non-perturbative gluon with energy fraction 
Xq, using a generating function 



f:/a^n,^^5(iogf -Elog 



After this, the azimuthal angles of each emission are selected according to a flat 
distribution and a number of correction factors are included to produce a weight 
for the generated partonic state. LDCMC can either produce weighted events or 
use the weight as a hit-or-miss probability to produce unweighted events. The 
weight can include many things such as: 

• The full splitting functions (also for quark propagators) instead of the sim- 
plified Optionally the full matrix element for a 2 — * 2 sub-collision 
can be used. In particular for the quark box closest to the virtual photon the 
full off-shell matrix element can be used. 

• The standard Sudakov form factors. 

• Emissions not satisfying the LDC constraint in eq.(]29|) are given zero weight. 

• The running of a^. The scale is taken to be the transverse momentum of 
the emitted parton which, due to the LDC constraint, is always close to the 
highest scale in the emission. 
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Name 


QCD cascade 


applicable 


processes 


event record 


c 

k_ 


48 




49 


>c 


forward evolution 
Willi ov^r ivi 


ep 


1*9* qq 


parton level 
weigiiLeu. eveiiLS 


( 


14 


3C 


]A 

54 


DE 


backward evolution 
with CCFM 

using unintegrated 
gluon density 


ep,'jp, pp 


1*9* ^ qq 
1*9* - QQ 

-yg* J/ipg 

9*9* qq 

n*n* (1(1 
9 9^ VV 


parton level 
unweighted events 
hadronization 
via Jetset/Pythia [ 132| 


L 

[ 


50 


1 


51 


1 


forward-backward 
symmetric LDC 
evolution 


ep 


* * — 

1 9 qq 
1*9* ^ QQ 


weighted or 
unweighted events 
final state cascade 
hadronization 
via Jetset/Pythia [|132| 



Table 4: Summary and overview over existing Monte Carlo event generators for 
small X physics. 



Averaging over weights for a given xq, x and Q^, it is then possible to fit the input 
gluon (and quark) distribution(s) to describe eg. F2(x, Q^). 

Using these input distributions, parton-level events can be generated. Final-state 
parton cascades are added using the colour-dipole cascade implemented in Ari- 
adne for e~^e~ annihilation but only allowing emissions which are below the LDC 
constraint in eg. (p9[) . Finally, standard Lund string fragmentation can be added 
to produce final-state hadrons. 

LDCMC is distributed together with the Ariadne program, available from 
http : / /www . thep . lu . se/~leif /ariadne but in using it, one should keep in mind, 
that the reproduction of data on eg. forward jets is very poor. The problem can be 
traced to the non-singular parts of the gluon splitting function. In a newer (not yet 
released) version, it is possible to allow only gluonic chains and only the singular 
parts of the gluon splitting function. The results are then consistent with what is 
obtained with Smallx and Cascade. 

6 Conclusions 

In this summary report we presented the state of the art of small x physics in 
the year 2001. Significant progress has been made in the understanding of the 
small X evolution equations, CCFM and BFKL. It has been possible for the first 
time to describe the structure fucntion F2{x,Q'^) and also hadronic final state 
measurements, like forward jet production, with the CCFM evolution equation 
implemented into a Monte Carlo program. In more detailed studies the need for 
improving the small x splitting functions became evident, also from considering 
next-to-leading corrections to the BFKL equation. However, considering these 
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improvements in detail, it became also clear, that this is not a trival task. 

Significant progress has also been made in the understanding of /c^-factorization 
in general and the calculation of the off-shell matrix elements. In certain ap- 
proximations the off-shell matrix elements are already calculated to order 0{al). 
However, some critical points still need to be clarified, as the gauge invariance of 
the /c^-factorization approach is not yet clear, if considered beyond leading order. 

For the first time, a comparison of all available paramerisations of unintegrated 
gluon distributions was made, showing significant differences, which indicate, that 
more (and more exclusive) measurements need to be used in constraining the un- 
integrated gluon distributions further. For the first time, we are in a position to 
try to perform a global fit, similar to those using the coUinear approximation, to 
determine the unintegrated gluon density. 

Our understanding of small-x physics is far from complete. There are a number 
of theoretical and phenomenological issues which need to be further settled. In 
this overview we have mentioned many such issues. On the theoretical side, we 
need to understand the convergence of the perturbative expansion of the evolution 
kernel, and we also need to calculate the impact factors to next-to-leading order. 
On the phenomenological side it is important to understand the importance of 
non-singular terms in the gluon splitting function especially when implemented in 
event generators. It is also important to understand the uncertainties involved with 
fitting unintegrated parton distributions to data. These parameterizations should 
then be compared to a wide variety of data, possibly requiring the calculation of 
new off-shell matrix elements for other than the mentioned processes. 

Finally, we want to note that not all aspects of small-x physics have been covered 
in this paper. Phenomena such as rapidity gaps, shadowing effects and multiple 
interactions are also interesting aspects of small-x evolution, and will be included 
in forthcoming meetings and publications. 

Appendix: The small- cc collaboration 

Today the work on these issues is spread out on a number of different small groups 
around the world. Although some of these group are already collaborating on an 
informal basis, it was agreed in the meeting, that the future work could be more 
coordinated. As a consequence of this a Small-x Collaboration was formed. The 
idea is to start small and informal and to set up a web site an a mailing list, but 
also to organize small meetings such as the one held in Lund resulting in this paper. 

The web site is located at http://www.thep.lu.se/Smallx and will among other 
things contain a compilation of subroutines implementing parameterizations of 
unintegrated parton distributions and off-shell matrix elements and compilations 
of relevant papers, theoretical as well as phenomenological and experimental. The 
mailing list would be used to announce new results, to ask for help from experts 
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in the field and so on. 



The Small-x Collaboration will, of course, be open to anyone in the field. To 
join the collaboration one can simply join the mailing list by sending a mail to 
smallx-subscribeOthep .lu.se (more instructions can be found on the web site). 
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